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ABSTRACT 

Recent observations in X-rays and gamma-rays of nearby FRI radio galaxies have 
raised the question of the origin of the emission detected in the termination structures 
of their jets. The study of these structures can give information on the conditions for 
particle acceleration and radiation at the front shocks. In addition, an evolutionary 
scenario can help to disentangle the origin of the detected X-ray emission in young 
FRI sources, like some Gigahertz Peaked Spectrum AGNs. This work focuses on the 
nature and detectability of the radiation seen from the termination regions of evolving 
FRI jets. We use the results of a relativistic, two-dimensional numerical simulation of 
the propagation of an FRI jet, coupled with a radiation model, to make predictions for 
the spectra and lightcurves of the thermal and non-thermal emission at different stages 
of the FRI evolution. Our results show that under moderate magnetic fields, the syn- 
chrotron radiation would be the dominant non-thermal channel, appearing extended 
in radio and more compact in X-rays, with relatively small flux variations with time. 
The shocked jet synchrotron emission would dominate the X-ray band, although the 
shocked ISM/ICM thermal component alone may be significant in old sources. Inverse 
Compton scattering of CMB photons could yield significant fluxes in the GeV and 
TeV bands, with a non-negligible X-ray contribution. The IC radiation would present 
a bigger angular size in X-rays and GeV than in TeV, with fluxes increasing with 
time. We conclude that the thermal and non-thermal broadband emission from the 
termination regions of FRI jets could be detectable for sources located up to distances 
of a few 100 Mpc. 

Key words: galaxies: jets-hydrodynamics-galaxies: evolution-X-rays: galaxies- 
gamma-rays: galaxies-radio continuum: galaxies 



1 INTRODUCTION 

Extragalactic jets from Active Galactic Nuclei (AGN) in- 
ject energy in the interstellar and intracluster media (ISM 
and ICM, respectively) at a rate between ~ 10 42 to 
~ 10 46 ergs _1 , depe nding on the source. F anaroff-Riley 
sources of type I (FRJ. lFanaroff fc RilevH"l974l ) fall typically 
on the lower edge of this power s pectrum. They show rela - 
tivistic velocities at parsec-scales (|Celotti fc Ghise llini 2008) 
and disrupted structure at kiloparsec scales, whereas the 
more powerful FRII jets keep collimated up to the medium 
interaction point, in which hot-spots can be observed at dif- 
ferent frequencies. The interaction of the jet with the am- 
bient in FRI and FRII galaxies could be important to the 
extent that AGN feedback has been claimed to be a possi- 
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ble solution for the cooling flow problem via shock-heating 
or mixing (e.g. Quilis, Bower & Balogh 200l|; IZanni et al.l 
120051 : iMcNamara fc NulserJl2007t ). Also, this interaction can 
give rise to heating and particle acceleration via shocks, in 
which thermal and non-thermal radiation is produced and 
can be used to study the p roperties of the flow and the 
mediu m. Following this idea, iHeinz. Reynolds fc Begelmanl 
(1 19981) used a simple evolutiona ry model based on the work 
bv lRevnolds fc Begelmanl l|l997l ) to obtain the X-ray bright- 
ness of the thermal emission for different initial jet proper- 
ties. They claimed that only for dense enough cluster me- 
dia, the count rates obtained would be enough to detect this 
emissi on even for powerful FRII iets. lKino. Kawakatu fc ltd 
i|2007l ) have also derived estimates for the thermal MeV 
emission from cocoons of radio galaxies depending on their 
age, with the result that only young cocoons, with ages 
<C 10 7 yrs could be detected at this energy band by present 
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space-observatories. IZanni et alj (|2003T l performed a series of 
simulations of supersonic and underdense jets in a decreas- 
ing pressure atmosphere and showed that jets evolve in two 
different phases regarding their high-energy thermal emis- 
sion: a phase in which the shell formed by shocked material 
is highly overpressured and radiative, and a later phase in 
which the shock is weaker and a deficit of X-ray emission is 
expected from the lobes. 

Regarding observations, iKataoka et al.l (|2003l) reported 
on the Chandra detection of faint, extended X-ray emis- 
sion &^m_th*jets_and_l£besc^he radio galaxy 3C 15 (see 
also lHarris fc Krawczvnskill2006h . This emission is spatially 
correlated with tha t observed at 8.3 GHz radio frequencies 
l|Leahv et al. 1 1 19971 ). The authors suggested that the same 
electron population responsible for the radio synchrotron 
emission upscatters t he CMB photons to produc e the dif- 
fuse X-ray radiation. ISiemiginowska et al.1 (|2008l ) reported 
the detection of X-ray emission from Gigahertz Peaked Spec- 
trum (GPS) and Compact Steep Spectrum (CSS) sources 
(13 quasars and 3 galaxies, all of them powerful sources), and 
claimed that this radiation is most likely related to the accre- 
tion power in all but one of the studied sources, in which the 
emission could be generated in the jet. They also discussed 
the possibility that the X-rays were produced in the bow 
shoc k formed by the exp andin g jet but found no evidence for 
this. iKraft et all (|2003h and ICroston. Kraft fc Hardcastlel 
(|2007 ) reported detection of X-ray emission in Cen A and 
NGC3801 using Chandra, which was interpreted in terms 
of the bow shock driven by the injection of a jet. Model- 
ing the emission as thermal, they obtained bow-shock Mach 
number s between 4 and 8. However, deeper o bservations of 
Cen A l|Kraft et al.ll2007l ; ICroston et al.ll2009r i showed that 
the emission from a bow-shock region around the south-west 
lobe is better interpreted as of synchrotron origin, implying 
that the shock is strong enough to accelerate particles up 
to Lorentz factors of ~ 10 8 . In other sources, like the radio 
galaxy Fornax A, the lobes seem to emit non-thermal X-rays 
through inverse Compton (I C) of cosmic microw ave back- 
ground (CMB) photons fe.g. iFeigelson et al.lll995T ). whereas 
the larg e-scale jet of M 87 woul d be also a synchrotron em it- 
ter (e.g. IWilsonll2003l ; see also lKataoka fe Stawarzl [20051 for 
a discussion on extended jet emission and possible origins). 
Recently, the detection by Fe rmi of extended Ge V emission 
in the radio lobes of Cen A l|Abdo et al.ll2~010al ). likely via 
IC scattering of CMB photons, shows that acceleration up 
to VHE is taking place in the disrupted jet region. 

iPerucho fc Martj (^OOTT l (PM07 hereafter) performed a 
simulation aim ed to test the F RI jet evolution parad igm 
i|BickneUlll984n and the model bv lLaing fc Bridle! |2002) for 
the FRI jet of the radiogalaxy 3C 31. The simulation was 
done using a numerical code for relativistic hydrodynam- 
ics, based on High-Resolution-Shock-Capturing schemes, to 
which it was added an equation of state that allows for a 
specific treatment of two families of particles, leptons and 
baryons, and computes the adiabatic index in terms of the 
composition of each cell. The jet was injected in the numer- 
ical grid at 500 pc from the active nucleus, with a radius 
of 60 pc. The ambient medium, composed by neutral hy- 
drogen, has a profile in pressure, density and temperature. 
Such a profile is required in most jet evolutionary models to 
account for the jet collimation at large distances. Further- 
more, all models in which the jet is decelerated shortly after 



being ejected, within distances 1-10 kpc from the nucleus, 
require a gradient in the ambient pressure that permits the 
jet to prevent disruption due to external mass loading. This 
profile includes the contribution from a core region, domi- 
nant for distances up to ~ 1.5 kpc, and a more extended, 
hotter and less dense contribution from the galaxy group, 
which dominates at large distances. For a detailed discus- 
sion on the X-ray properties used to cha racterize the ex- 
ternal med ium in FR-I sources see, e.g., lHardcastle et all 
|2002l) and lLaing fc Bridle! (|2002l ). 

The jet, leptonic in composition, was injected with 
a velocity Wjo = 0.87 c, density ratio with the ambient 
Pjo/PaO = 10 , pressure ratio with the ambient Pjo/P a o — 
8, and temperature 4 x 10 9 K, resulting in a kinetic lumi- 
nosity Lj = 10 44 erg s _1 . The simulated jet evolved during 
~ 7 x 10 6 yrs up to a distance of 15 kpc. For further details 
on this simulation, we refer the reader to PM07. 

In the present work, we use the results from this simula- 
tion of an FRI jet interacting with the ISM/ICM to compute 
the produced thermal and non-thermal emission for different 
source ages. We have coupled a simplified radiation model 
for the cocoon an d the shell applied already to the context 
of microquasars (|Bordas et al.1 12009| ) to the results of the 
simulations of PM07. In this work, we make a specific use of 
the terms cocoon and shell. Namely, we refer to the cocoon 
as the region of the jet shocked material, starting already 
at the recollimation (see below), and the shell as the region 
of shocked external medium. We have also covered source 
ages older than 7 x 10 6 yrs using extrapolations of the main 
hydrodynamical paratemers derived from the simulation re- 
sults in PM07. In this way, we can make predictions for the 
flux and the spectral evolution of the thermal (X-rays) and 
non-thermal (radio to gamma-rays) emission of an FRI jet 
for a broad age range: 10 5 — 10 8 yr. We discuss the relevance 
of the thermal and the non-thermal radiation, and the pos- 
sibility to produce high- (HE) and very high-energy (VHE) 
from the termination regions of FRI jets. 

The paper has been organized as follows: in section [5] we 
present the emission model and its results for different stages 
of the cocoon and shell evolution, characterized using the 
simulations of PM07. The discussion of the results and the 
conclusions are presented in Sect.[3]and Sect.|4j respectively. 



2 THERMAL/NON-THERMAL EMISSION 
FROM FRI JETS 

2.1 The model 

The model adopted here to study the non-thermal emis- 
sion of the termina t ion si te of an FRI jet has been adapted 
from iBordas et ail (2009), in which the non-thermal radi- 
ation of a microquasar jet termination region was studied. 
In that paper, the dynam ics w as ba s ed on the works by 
iKaiser fc Alexander] l|l997ft and|Falle] (|l99ll ). whereas here 
the dynamics has been extracted from the simulations by 
PM07. The thermal emission has been computed using the 
information on the density and temperature obtained from 
these simulations. For details on the properties of the jet at 
injection and the external medium, we address to Tables 1 
and 2 in PM07 (see also Table 3 of the same work for a com- 
parison with Cen A and the radiogalaxy NGC 3801). Since 
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in the present case the jet is disru pted and a strong re verse 
shock is not produced (unlike in iBordas et al.l |2009| ). the 
shock in the jet (reverse shock) has not been considered. In- 
stead, we have accounted for the strong recollimation shock 
as the particle accelerator in the cocoon. For illustrative pur- 
poses, we present in Fig. [1] a density map, with isobars, of 
the simulated jet after 7 x 10 6 yr of evolution (see PM07). 
The inset in Fig. [T] shows a zoomed view of the head of the 
jet. 

2.1.1 Non-thermal particles and their emission 

Most of the accelerated electrons do not cool down sig- 
nificantly before overcoming the jet disruption point (DP) 
and reach later on the turbulent cocoon. This is due to the 
fast motion downstream the recollimation shock, which ef- 
ficiently carries particles away along the jet axis down to 
the DP. Also, the compression of the shocked jet walls by 
the cocoon material prevents expansion and therefore adi- 
abatic cooling until electrons reach the DP. Beyond that 
point, relativistic electrons spread in the cocoon via turbu- 
lent advection and diffusion. 

Since the material becomes trans- or subsonic after the 
recollimation shock, the pressure in the shocked jet, the co- 
coon and the shell should be relatively similar, the shell be- 
ing denser but colder than the cocoon. In addition, the dom- 
inant photon field is the (homogeneous) CMB one, with ra- 
diation energy density mcmb » 4.2 x 10" 13 (1 + z) 4 erg cm -3 
(we use z = in our calculations), provided that the emitter 
is located at a distance > 2.6 (-Lnuc/10 43 erg s _1 ) kpc from 
the galaxy nucleus (L nuc is the nucleus luminosity). This al- 
lows us, at this stage, to simplify the cocoon region as an 
emitter with homogeneous properties (one zone) with the 
recollimation shock as the injector of accelerated particles. 
The external medium shocked by the bow shock, i.e. the 
shell, should be also mostly subsonic, and the same consid- 
erations regarding the photon field apply there. Therefore, 
we have also adopted an homogeneous emitter approxima- 
tion for the shell. 

In both the shell and the cocoon, the luminosity in- 
jected in the form of non-thermal particles has been taken as 
a 10% of the total jet kinetic luminosity, i.e. L nt = 0.1 Lj = 
10 43 erg s _1 . The magnetic field B has been fixed taking 
the magnetic energy density, ub = B 2 /8iv, to be 10% of 
the ram/thermal pressure. Concerning particle acceleration, 
the recollimation shock has been assumed to be relativis- 
tic, with an acceleration rate E = r/qBc with r\ = 0.1. 
For the bow shock, we have adopted the expression for a 
non-relativistic shock, in which r\ = ^(ub/c) 2 (e.g. Drury 
1983), where Vbs is the bow-shock velocity (typically here 
of ~ (1 — 2) x 10 8 cm s _1 ). These acceleration rates are 
to be compared to the syn chrotron and IC loss rates (e.g. 
iBlumenthal fc G ould 1970]) to derive the maximum energy 
of electrons. The cooling timescales (t C ooi = —E/E coo i) of 
synchrotron and (Thomson) IC processes are: 

t syn r;4x 10 12 {B/WfxG)~ 2 (S/lTeV)" 1 s , (1) 

and 

tic « 1.6 x 10 13 (u/KT 12 erg cm^ 3 )" 1 (E/lTeVy 1 s , (2) 

respectively, being u the total radiation energy density. An 
escape time has also to been considered since particles with 



enough energy would escape the accelerator. This is derived 
by taking the gyroradius of the most energetic particles 
equal to the size of the accelerator, i.e. the recollimation 
and bow-shock widths (Hillas 1984). We do not consider the 
possible role of Fermi II stochastic or shear a cceleration in 
the d isrupted jet and cocoon regions (see, e.g.. iRieger et al.l 
l2007t see also lO' Sullivan et alj|2009l for a deeper analysis of 
Fermi II particle acceleration in the context of the Lobes of 
Cen A), although these processes may be absorbed by our 
phenomenological treatment of the particle acceleration in 
the cocoon. 

The properties of the non-thermal emitters in the co- 
coon and the shell are characterized by the ram/thermal 
pressure and the bow-shock velocity (and the shock sizes 
when cooling is unefficient), which determine the magnetic 
field, the synchrotron emission, indirectly the IC emission, 
and the acceleration efficiency. These conditions have been 
parameterized making use of the results of the simulations 
of PM07 and their extrapolation to earlier and later times, 
covering an age range t BTC = 10 s — 10 s yr. We do not expect 
significant uncertainties from the extrapolations as long as 
the medium properties present the same properties at larger 
distances than those covered by the simulated jet. Some of 
the model parameters are listed in Table[T]for both the shell 
and the cocoon regions. 

The spectral aging of the non-thermal particle popu- 
lations has been modeled considering the evolution of the 
physical conditions in each interaction region. The parti- 
cle energy distribution at a given time, N(E,t SIC ), is calcu- 
lated by adding the different evolved injected populations, 
Q(E,t) (oc E~ p ), from t — up to t BTC . The time resolu- 
tion of particle injection is At < t coo i(t). Maximum parti- 
cle energies, -Emax(t), are also computed for each time step 
due to the time dependence of the magnetic field, the ac- 
celerator size and the shock velocity. For simplicity, a spec- 
tral index p = 2.1 has been used in our calculations for 
both the recollimation and the bow shock. We note that, 
together with synchrotron and IC cooling, the expansion 
of the jet termin ation structure intr oduces an adiabatic 
loss timescale (see IBordas et al. 2009) ~ l/v h « (5/3) two, 
where I is the size of the whole structure. In addition to 
synchrotron and IC process es, relativistic Bremsstrahlung 
ijBlumenthal fc Gouldl ll970T ) could also take place in the 
shell, and protons may be accelerated and eventually could 
interact with the s hocked jet medium t hrough proton-proton 
(pp) collisions fsee lKelner et al1l2006l ). However, the densi- 
ties n of targets for relativistic Bremsstrahlung and pp emis- 
sion in the shell, the largest in the jet termination region, 
are low, and the cooling timescales: 

trel.br/pp ~ 10 18 (n/10 -3 cm" 3 ) s > t 8rc . (3) 

This implies radiation efficiencies much smaller than 
those of synchrotron and IC. In additio n, it is worth not - 
ing that synchrotron proton emission (e.g. Aharonian 2000), 
under equipartition magnetic fields, could overcome IC radi- 
ation around 100 MeV, although only for very young sources 
(^C 10 5 yr) this component may be significant. We do not 
further consider either relativistic Bremsstrahlung or proton 
radiation processes in this work. 



4 P. Bordas, V. Bosch-Ramon and M. Perucho, 




3.0X10 3 6.0X10 3 9.0X10 3 1.2X10 4 1.5X10° 1.8X10 4 



r(pc) 

Figure 1. Snapshot of the density in color scale for the simulation of an FRI jet after 7 X 10 6 yr (for details, see PM07). The shell is 
clearly seen, as well as the recollimation shock, the disruption point and the turbulent cocoon region. The inset shows a zoom around the 
head of the bow shock. Pressure contours at the level of 2 and 8 xlO -12 erg cm -3 are labeled. The highest pressure zones correspond to 
the regions close to the Rankinc-Hugoniot conditions and are used to compute the RH thermal emission, whereas the Av contribution 
makes use of the values averaged over the whole shell (see text for details) 



Table 1. Model parameters for the shell and cocoon for three different source ages used to compute the thermal and non-thermal 
emission. 



Parameter 


t S rc = 10 5 yr 


t src = 3 x 10 6 yr 


t src = 10 s yr 


Bow shock velocity vi, (c) 


8.7 X icr 3 


6.2 x 10~ 3 


4.4 x 10~ 3 


Shell density p sri (g cm -3 ) 


1.0 x icr 24 


4.0 x 10~ 26 


3.4 X 10~ 27 


Shell temperature T s h (K) 


8.3 x 10 7 


2.7 x 10 7 


3.5 x 10 6 


Shell radius r s i, (cm) 


1.6 x 10 20 


6.5 x 10 21 


1.2 x 10 23 


Shell and cocoon magnetic field B (G) 


2.5 x 10~ 4 


3.1 x 10~ 5 


6.4 x 10~ 6 


Cocoon radius r coc (cm) 


5.5 x 10 20 


2.2 x 10 20 


4.0 x 10 20 


Shell maximum energy -E^ ax (TeV) 


17.0 


33.5 


51.0 


Cocoon maximum energy S^ax (TeV) 


1.5 x 10 3 


4.4 x 10 3 


9.1 x 10 3 



2.1.2 Thermal emission 



The thermal emission has been computed making use of 
the simulation results and their extrapolation to the i S rc- 
range considered here. Given the strong density dependence 
of thermal Bremsstrahlung, we have only accounted for 
the contribution from the shell, much denser than the co- 
coon. Furthermore, we have simplified the calculations of the 
thermal radiation as it would be coming from two regions 
(see Fig. [TJ. One, cooler (ultraviolet -UV-/soft X-rays) but 
brighter, corresponds to the averaged shell conditions {Av), 
and another one, fainter but hotter (hard X-rays), corre- 
sponds to a region close to the apex of the bow shock (RH), 
in which the shell has properties close to those given by the 
jump conditions of Rankine-Hugoniot. The volume of the 
latter region is about 3 — 4% of that of the whole shell (see 
inset in Fig. [JJ, which corresponds to the volume limited 
by the isobars satisfying P > | Pr-h, where Pr-h is the 
shell pressure right behind the bow shock and corresponds 
to the Rankine-Hugoniot jump conditions. Given the high 
temperatures in the shell, we have calculated the thermal 
Bremsstrahlung assuming that the plasma is fully ionized, 
with electrons and protons in equipartition. At this stage, 
we have not considered line emission. 



2.2 Results 

We have studied the thermal and non-thermal emission pro- 
duced in the shell and the cocoon separately. We consider 
here the contribution from a single FRI jet, so the predicted 
luminosities should be scaled by a factor of two to obtain 
the whole source emission under similar ambient conditions 
for jet and counter-jet. 

2.2.1 Non-thermal emission 

The bow shock can accelerate electrons up to energies 
£max ~ 17 TeV (10 5 yr) to 51 TeV at (10 8 yr), and is 
limited by synchrotron losses at any time. This trend of 
higher E m&K for older sources comes from the energy gain 
to loss ratio oc B~ l under synchrotron dominance. Since Vb 
decreases moderately, from « 2.8 x 10 s to 1.3 x 10 s cm s _1 , 
the strongest variation in the shell acceleration rate comes 
from the P-evolution, which goes from s=s 2.5 x 10~ 4 (10 5 yr) 
to 6.4 x 10" B G (10 8 yr). The shell IC emission is dominated 
by scatterings with CMB photons. In the cocoon region, we 
have assumed the recollimation shock to be the accelera- 
tor site. The maximum energy also grows here, going from 
« 1.5 x 10 3 to 9 x 10 3 TeV. Since Mb is proportional to the 
pressure and the latter is similar in the cocoon and the shell 
(see Figs. 5 and 6 in PM07), B is also similar in both re- 
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Figure 2. Computed SEDs of the non-thermal synchrotron and IC emission from the shell and the cocoon for three different ages: 10 5 
(solid line), 3 X 10 6 (long-dashed line) and 10 s yr (dotted line). We have accounted for the contribution of one jet/medium interaction 
region only, so the values here should be scaled by a factor of two to get the whole source emission. 



gions. As in the shell, synchrotron losses dominate for the 
magnetic field strengths and ages considered here. The high 
values of E max in the recollimation shock are expected since 
the acceleration rate is assumed to be ~ (c/v s ) 2 times more 
efficient here than in the non-relativistic bow shock. The 
large distance of the recollimation shock to the galaxy nu- 
cleus makes the CMB IC to dominate over other IC compo- 
nents, although for very young sources the galaxy nucleus 
could be relevant. 

The non-thermal spectral energy distributions (SED) 
for the cocoon and the shell, at t Brc — 10 s , 3 x 10 6 and 
10 s yr, are shown in Fig. [2] The obtained radio and X- 
ray synchrotron luminosities in both regions are at the level 
of 2 x 10 41 erg s _1 . The approximate constancy of the lu- 
minosities with time is due to the fact that particles have 
reached the steady state at t src through synchrotron coolinj£] 
and the assumed constancy of L nt . The decrease of B with 
time, and therefore the growth of t syn , is compensated by 
the increase of time available for cooling. The synchrotron 
break frequency, corresponding to the electron energy at 
which t syn (E) « t SIC , and the highest synchrotron frequency, 
Vsyn max oc -B-E max , are shifted down for older sources. The 
former effect makes the radio luminosity to increase at the 
late stages of the evolution of both the cocoon and the shell, 
whereas the latter decreases the X-ray luminosity in the shell 
due to the decrease of u syn max with time. The slightly dif- 
ferent conditions in the shell yield a higher break frequency, 
which implies a factor ~ 2 lower radio luminosity in this 
region compared to that of the cocoon. 



1 Actually, the adiabatic cooling, as approximated here, takes 
~ 1/2 of the particle energy after a time ~ t src , the rest of the 
energy going to radiation. 



The IC luminosity grows as long as this process be- 
comes more efficient compared to synchrotron and adiabatic 
cooling, which is shown by the decrease of us/u ra d from 
R! 5 x 10 3 (10 s yr) to 4 (10 s yr). As expected from the 
Smax-values given above and the similar energy budget, the 
cocoon and the shell have similar HE luminosities, but the 
cocoon is few times brighter at VHE than the shell due to its 
much higher maximum frequency. In both regions the bolo- 
metric IC luminosities grow similarly with time, reaching 
~ 10 42 and 10 41 erg s" 1 at HE and VHE, respectively. 

The lightcurves for the luminosities in radio (5 GHz x 
Ls GHz), X-rays (bolometric: 1-10 keV), HE (bolometric: 
0.1-100 GeV) and VHE (bolometric: 0.1-100 TeV), for both 
the cocoon and the shell, are presented in Fig. [3] The 
lightcurves show in more detail the time behavior of the non- 
thermal radiation at different wavelengths discussed above. 
The complex and smooth shape of the lightcurves, most 
clear for the HE and the VHE emission, is a consequence of 
the complex hydrodynamical evolution of the whole inter- 
action structure propagating in an inhomogeneous external 
medium. 

2.2.2 Thermal emission 

As mentioned in Sect. l2.L2l thermal Bremsstrahlung is also 
expected from the shell. Figure U shows the three SEDs 
(t src = 10 5 , 3 x 10 6 and 10 8 yr) computed adopting a simpli- 
fied model for the thermal emitter of the shell considering 
the shell averaged values, Av, and the conditions right be- 
hind the bow shock apex, RH. The slowdown of the bow 
shock and the velocity dependence of the postshock temper- 
ature, oc Wb a , leads to a decrease in the peak of the thermal 
emission with time, whereas the increase of the shell mass 
yields higher thermal bolometric luminosities as the source 
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gets older. For the age range t BIC = 10 5 — 10 8 yr, the thermal 
luminosities go from 10 39 to few times 10 41 erg s _1 , with the 
shell and the hot postshock region components peaking from 
soft X-rays to UV and from hard to soft X-rays, respectively. 
The hot postshock region dominates the SED in hard X-rays 
by a factor of a few over the shell thermal and both shell 
and cocoon non-thermal components for t BTC = 10 8 yr. 



Figure [5] shows the thermal lightcurves (bolometric) for 
tsrc = 10 5 — 10 8 yr. Thermal Bremsstrahlung increases from 
tsrc = 10 5 to ~ 10 6 yr, when it shows a relative maximum. 
Then the luminosity slightly decreases until t BTC ~ 3 x 10 6 yr, 
time in which there is a transition in the external medium, 
from the denser galaxy core to the rarefied galaxy group 
medium (see PM07 for details). Later, the emission increases 
again. The component Av dominates the thermal bolomet- 
ric luminosity in young sources, but the component RH be- 
comes similarly bright at t BTC ~ 10 8 yr. 



3 DISCUSSION 
3.1 Radio 

The cocoon and the bow shock show both a similar pat- 
tern of their non-thermal radio emission, although a higher 
£ max makes the cocoon emission to extend to higher en- 
ergies. The accumulation and aging of the injected parti- 
cles lead to, for t Brc = 10 s yr, a break in the synchrotron 
spectrum around the radio frequencies. The cocoon would 
be the dominant radio emitter, with fluxes as high as ~ 
lCT 12 (d/100 Mpc)" 2 erg cm" 2 s" 1 or ~ 10 Jy at 5 GHz 
from a region of few times 10' (d/100 Mpc) -1 angular size. 
The spectral index would appear inverted due to particle 
aging, with a ~ 1 (F v oc v~ a ). The properties of the 
radio emission from the interaction jet-medium structure 
are comparable with those observed for instance in 3C 31, 
with radio luminosities at 4.75 GHz of about 3 x 10 40 erg/s 
l|Andernach et aflll992h . or with the ones of 3c 15, in which 
fluxes of a few xlO" 14 erg cm" 2 s" 1 (d ~ 300 Mpc) are 
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Figure 4. Computed SEDs of the shell thermal emission for three 
different ages: 10 s (solid line), 3 X 10 6 (long-dashed line) and 
10 8 yr (dotted line). The two components are shown, one corre- 
sponding to the whole shell (thick lines) , and another one related 
to a shell region with conditions similar to those of Rankine- 
Hugoniot (thin lines). 

found (jKataoka et all 120031 ) . The predicted shell radio flux 
is sligthly below the cocoon one, although effects of limb 
brightening may enhance the detectability of the former. 
The radio lightcurve is quite steady, with some small vari- 
ations. Despite the fact that the magnetic field gets weaker 
with time, the accumulation of radio emitting electrons com- 
pensates it, and the final emission keeps roughly constant. 

3.2 X-rays 

Thermal X-rays are produced in the shell, with a tem- 
poral evolution smoother and more complex than in the 
case of a homogeneous medium. The existence of differ- 
ent emitting regions in the shell would lead to a rela- 
tively flat thermal spectrum in X-rays, with a bolometric 
flux from ~ 10" 15 (d/100 Mpc) -2 (10 5 yr) to few times 
1CT 13 (d/100 Mpc)" 2 erg cm" 2 s" 1 (10 8 yr). We note that 
the thermal emission would be restricted to different an- 
gular size regions depending on the photon energy. The 
hard X-ray photons would come from the apex of the bow 
shock, with typical angular size of a few 1' (d/100 Mpc) -1 
(for t STC ~ 10 8 yr), and the lower energies would be dom- 
inated by the whole cooler shell emission, with an angu- 
lar size of a few 10' (d/100 Mpc) -1 . Limb brightening ef- 
fects could play a role, showing a thin structure along the 
limb of the shell with the hottest region at the apex. It 
is worth noting that under the adopted L nt -value and B- 
equipartition fraction, the shell thermal emission dominates 
the emission except in hard X-rays. We also note that 
given the moderate velocities of the bow shock, thermal 
p hotons cannot reach energies as high as those discussed 
in lKino. Kawakatu fc ltd (|2007l ). Nevertheless, for the shell 




le+38 1 ' — ' — ' — 

le+05 le+06 le+07 le+08 

Figure 5. Evolution of the computed thermal bolometric lumi- 
nosity with time. Two components are shown, one computed with 
density and temperature averaged in the whole shell, Av (dotted 
line), and the other one corresponding to a region with condi- 
tions similar to those of a strong shock, RH (long-dashed line). 
The summation of both components is also shown (thick solid 
line). 

properties considered in this work, the thermal cooling time- 
scale tth ~ 2.5xl0 9 ( T ^) a5 ( 10 _? c - m _, ) -1 yr is greater than 
t SIC . Assuming that the bow shock keeps being adiabatic and 
strong all along the source age, and hence not displaying a 
transition to a much weaker shock regime (see e.g. Zanni et 
al. 2003), the thermal bolometric luminosity increases with 
time. 

Regarding non-thermal X-rays, the dominant emis- 
sion comes also from the cocoon, with fluxes ~ 
10" 13 (d/100 Mpc) -2 erg cm" 2 s~\ although again limb 
brightening effects may increase the shell detectability. In 
fact, in the case of Cen A, the shell seems to be the dom- 
inant source of non-thermal X-rays (Croston et al. 2009). 
This difference could be explained by a higher E max in the 
shell of that source. In addition, a relatively recent de- 
crease in jet power would have affected first the cocoon 
synchrotron emission, making this radiation fainter while 
the shell emission would remain at similar levels for a time 
> lOkpc/c ~ 3 x 10 4 yr. The lifetime of X -ray synchrotron 
electrons, ~ 10 11 s, is much shorter than in radio, and <C t arc 
as well. This implies that these particles may not have time 
to reach the whole emitting structures, and their radiation 
may come mostly from the inner regions of the cocoon or 
the bow-shock apex. Note that this may lead to a violation 
of the assumption of an homogeneous emitter. This X-ray 
synchrotron emission concentrated around the recollimation 
shock is compatibl e with the large-scale je t X-ray emission 
found in 3C 31 bv lHardcastle et al.l l|2002h . If a strong rec- 
ollimation shock is indeed the origin of these large-scale jet 
X-rays, then the hypothesis that jet disruption in 3C 31 
is caused by shock triggered instabiliti es is favored against 
stellar wind mass-load (as proposed bv lLaing fc Bridleil2002l ; 
see also PM07). 

Like the radio emission from the cocoon, particle aging 
makes the X-ray synchrotron spectrum in the cocoon and 
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the shell to be relatively soft, with photon index F > 2 (n v oc 
v~ T ). The total non-thermal X-ray flux is roughly constant 
for the explored range of t S rc , although the shell contribution 
decreases significantly with time. Close to the recollimation 
shock, before reaching the DP, the spectrum in soft X-rays 
would appear harder, since the corresponding emitting elec- 
trons could not have time to cool. We note that X-ray fluxes 
at ages ~ 10 6 — 10 7 yr wo uld not be very different from those 
found in young sources bv lSiemiginowska et al.l (|2008f ) . Also, 
the non-the rmal X-ray fluxes ~ 10 ~ 14 erg cm -2 s" 1 of 3C 15 
reported bv lKataoka et alj l|2003l ) (see Fig. 8 in that work) 
imply a non-thermal luminosity of ~ 1.3 x 10 41 erg s" 1 at 
300 Mpc, in good agreement with the values predicted here. 

As noted in Sect.[T] the extended emission in some FRI 
jets has been interpreted as IC instead o f synchrotro n emis - 
sion. As shown for instance in Fig. 5 of llsobe et al. I (120061 ) 
for Fornax A, in some cases the X-ray emission is difficult to 
reconcile with a synchrotron origin. We note however that 
the predicted IC X-ray s pectra (Fig, wo uld appear sim- 
ilar to the one shown in llsobe et al. (2006). Therefore, for 
sources older than those considered here and/or lower B- 
values, in which synchrotron emission would be less rele- 
vant, cocoon and shell IC would dominate the X-ray output 
(covering roughly the whole cocoon/shell structure given the 
long cooling timescales of IC X-ray electrons). We remind 
that the complex medium of the turbulent cocoon region, 
not treated here, may enhance the X-rays in certain com- 
pact regions. Much lower surface brightness could prevent 
the detection of the more diffuse X-rays coming from larger 
regions of the cocoon. X-ray emitting electrons may also lose 
a significant fraction of their energy right after the recolli- 
mation shock. Given the transrelativistic velocities in that 
region, Doppler boosting would beam the emission favoring 
certain viewing angles. This effect has not been accounted 
for in the present study. 

Finally, we remark that thermal emission seems un- 
avoidable given the medium densities and bow-shock tem- 
peratures, whereas IC depends on L nt , as well as syn- 
chrotron, which also depends on B, none of these two pa- 
rameters being well determined. 



3.3 Gamma-rays 

The predicted SEDs in the HE-VHE range are similar for 
both the cocoon and the shell, although the latter shows 
a lower maximum photon energy. The HE SED is close to 
flat, and becomes steeper at VHE. We have not accounted 
for EBL gamma-ray absorption, which would become sig- 
nificant at distances larger than 100 Mpc. Regarding the 
lightcurve, the emission increases with time mainly due to 
the increasing efficiency of the CMB IC channel as com- 
pared with synchrotron and adiabatic losses. We note that 
the gamma-ray fluxes for a source with i arc ~ 10 8 yr are 
around ~ 10~ 12 (d/lOOMpc) erg cm" 2 s _1 . At HE, such a 
source may require very long exposures to be detected by, 
e.g., Fermi, although it cannot be discarded that very nearby 
sources, or sources with bigger non-thermal efficiencies or jet 
powers, could be detected after few years of observations. 
Actually, Fermi has already detect ed several FRI gal axies 
up to few hundred Mpc distances (|Abdo et al.ll2010bh , in- 
cluding t he extended radio lobes of Cen A (at a distance 
~ 4Mpc; |lsrae]||l998h . 

presenting fluxes similar to those pre- 



dicted here. At VHE, the fluxes would be detectable by the 
current instruments, although the extension of the source, 
of tens of arcminute at 100 Mpc, and the steepness of the 
spectrum above ~ 100 GeV, may make a d etection difficult. 
In th e case of Cen A, detected by HESS l|Aharonian et all 
2009), the emission seems to come only from the core, but 
this is expected given the large angular size of the lobes of 
this source, which would dilute its surface brightness too 
much. In general, long exposures with present Cherenkov 
instruments, like HESS, MAGIC II and VERITAS, and the 
forthcoming CTA, may allow the detection of VHE emission 
from FRI jet lobes, and possibly carry out energy-dependent 
morphological studies. 



4 CONCLUSIONS 

We have applied a radiative model to a prototypical FRI 
jet characterizing the flow with the results of hydrodynami- 
cal simulations. Thermal Bremsstrahlung (X-rays), and non- 
thermal synchrotron (radio- X-rays) and CMB IC (X-rays- 
gamma-rays) have been considered as the emission mecha- 
nisms. 

From our study we conclude that, for moderate non- 
thermal luminosities, radio lobes of FRI radio galaxies are 
good candidates to be detected in the whole spectral range, 
with the radiation appearing extended in most of the energy 
bands. The precise extension of the emitting regions is hard 
to calculate, and may depend e.g. on the source distance, 
the instrument resolution and the capability to disentangle 
the non-thermal emission from the background contribution. 
Our study does not aim to provide specific values of the 
source extension at different energies, and offers only rough 
estimates of the overall emission morphology. 

Our results show that soft X-rays may be likely dom- 
inated by synchrotron emission up to ages ~ 10 8 yr, with 
IC tending to be dominant for older sources. Thermal X- 
rays seem unavoidable and may dominate in hard X-rays in 
old sources even if a non-thermal component is present. The 
low surface brightness may require long observation times 
for the detection in X- and gamma-rays, although the steady 
nature of these sources can help in this regard. Moderate res- 
olution radio and X-ray data, with long enough exposures, 
can allow the direct comparison between predictions of sim- 
ulations and observational data, thus giving clues on the 
hydrodynamics of the present flows and their surroundings. 
Also, any nearby galaxy of this kind can be a suitable can- 
didate for an eventual gamma-ray detection. Non-thermal 
synchrotron X-rays and HE and VHE gamma-rays provide 
suitable information to study particle acceleration in the jet 
termination regions. 
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